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Abstract 

We formulate and analyse a saturation model for the total 77 and 7*7* cross-sections 
and for the real photon structure function F2 {xyQ"^). The model is based on a picture 
in which the 7*7* total cross-section for arbitrary photon virtualities is driven by the 
interaction of colour dipoles, into which the virtual photons fluctuate. The cross-section 
describing this interaction is assumed to satisfy the saturation property with the saturation 
radius taken from the Golec-Biernat and Wiisthoff analysis of the j*p interaction at 
HERA. The model is supplemented by the QPM and non-pomeron reggeon contributions. 
The model gives a very good description of the data on the 77 total cross-section, on the 
photon structure function F2{x, Q"^) at low x and on the 7*7* cross-section extracted from 
LEP double tagged events. Production of heavy flavours in 77 collisions is also studied. 
Predictions of the model for the very high energy range which will be probed at future 
linear colliders are given. 



1 Introduction 



The concept of parton density saturation in high energy scattering mediated by strong in- 
teractions is certainly one of the most interesting recent developments in QCD theory and 
phenomenology p|-p!6|. It is well known, that the cross-sections for processes characterised by 
a large scale exhibit a steep power-like rise with the collision centre of mass energy squared 
W'^ , when W"^ ^ Q^. Such a rise, attributed to the so called hard pomeron [0,|T3) if continued 



to arbitrarily large energies would eventually lead to break-down of the S'-matrix unitarity. To 
avoid the apparent unitarity violation, effects related to screening (shadowing) phenomena have 
to be considered. They correspond to multiple bare pomeron exchanges and multi-pomeron 
interactions, and tame the steep rise at high energies. In the language of parton densities, 
the cross-section rises as the number of partons grows in the target, due to gluon emissions 
into the available rapidity interval, Y ~ InCW"^ / Q"^) . This growth may be described as evolu- 
tion in rapidity, with the evolution length Y. However, apart from creation of new partons a 
competing phenomenon of gluon recombination occurs, which reduces the number of partons. 
The recombination becomes increasingly important at high parton densities, i.e. for a large evo- 
lution length Y. This qualitative picture has a solid theoretical basis, well rooted in QCD [p]-[T5[|. 



A very non-trivial feature which arises from those studies is that the characteristic rapidity 
evolution length Y, for which the unitarity corrections become important, depends on the hard 
scale Q^. This statement may be inverted, leading to a notion of the y-dependent saturation 
scale Qs{Y), the characteristic scale for the transition between colour transparency and satu- 
rated cross-section regimes, at given rapidity Y. This phenomenon has recently been thoroughly 
studied through the measurement of ep inelastic scattering at HERA within a broad range of 

varying from the DIS large region down to the real photoproduction limit ~ ||18|| . 
Here = —q"^, where q is the four-momentum transfer in the process ep e'X. Since this 
process is largely controlled by one photon exchange, corresponds to the photon virtuality. 
Measurement of ep inelastic scattering permits determination of the virtual photon - proton 
total cross-section cr^*p{Q'^, W"^) for all virtual photon polarisation states, with W being the 
photon - proton collision energy. 

Golec-Biernat and Wiisthoff (GBW) managed to fit these data in a model incorporating 
the saturation property, with rapidity dependent saturation scale |T6|. In this case the 7*p 
cross-section is described in terms of the qq colour dipoles which the (virtual) photon fluctuates 
to according to a known wave function. The dipoles scatter off the proton with a cross-section 
which exhibits the colour transparency and the saturation property in the limit of a small dipole 
size, r <^ l/Qs(Y), and a large dipole size, r ^ l/Qs{Y), respectively. It is an encouraging 
result that the predictions of the simple model agree well with all the large rapidity data rang- 
ing from photoproduction to large Q^. This model would therefore provide a description of the 
transition between the soft and hard high energy scattering in QCD. The same model explains 
also properties of the cross-section for hard diffraction at HERA, in particular the ratio of the 



total to diffractive cross-section being constant with energy [19 



Thus, it is important to perform other tests of the GBW model which probe its universality. 
Two virtual photon interactions at high energies offer an ideal opportunity for such studies since 
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virtualities of both photons can vary, so that the properties of the model may be studied more 
extensively. Beside that, in the GBW model, the photon wave function is known, contrary to 
the proton wave function, so the two-photon data may be used to constrain the dipole-dipole 
cross-section itself. 



There exist several models of two-photon interactions which aim at describing the variation 
of the dynamics depending on the photon virtualities [^0|-|28|1. Most of those models combine 
the Vector Meson Dominance with the Parton Model suitably extended to the region of low 
virtualities. They do also usually rely on the Regge pole description of the high energy be- 
haviour of the total cross-sections. Some of the models describing the total cross-sections of 
real photons explore the minijet production mechanism [^,0. The 7*7* interactions have 
also been described within the dipole picture P5|-p5|, but possible saturation properties of the 
dipole-dipole cross-section have not been studied so far. The saturation ansatz may be useful 
to better understand the features of the available two-photon data and to formulate some in- 
teresting predictions for the two-photon physics at future linear colliders. 



In this paper we construct a generalisation of the saturation model for the two-photon case, 
compare its predictions with the experimental data and discuss the implications. The content 
of our paper is the following. In the next section we recall the GBW saturation model for 7*p 
scattering and present its formulation for 7*7* high energy interactions. We point out that by 
a suitable choice of the quark masses, the model can be used to describe the total 77 cross- 
section for two real photons, the 7*7* total cross-section for two virtual photons measured in 
double tagged e^e^ e^e~ +X events and the photon structure function F2 (x, Q^) of real (or 
virtual) photon for low values of the Bjorken parameter x. We show that in the region where 
the saturation effects become important the model gives steeper dependence of the cross-section 
on the collision energy than that obtained in the case of 7*p scattering. In Section 3 we present 
comparison with the available experimental data on a-y-y, cr^.^*, at low x and for 

the cross-section describing heavy flavour production in 77 collisions. Section 4 contains our 
predictions for the above quantities in the very high energy regime, which can be available in 
future linear ee, je or 77 colliders. Finally in Section 5 we present a summary of our results. 



2 The saturation model 

2.1 The Golec-Biernat— WiisthofF model 

The study of the total virtual photon - proton cross-section in the high W limit and for 
ranging from small to large values allows probing the transition from large to short distance 
physics in high energy scattering. Numerous analyses exist which study this transition |[18| , |3T[| . 
Among them, a very successful description is provided by the saturation model developed by 
Golec-Biernat and Wiisthoff |T^, in which the 7*p scattering is viewed upon as the scattering 
between {qq)dipoie and the proton. The colour dipoles (gg)dipoie represent virtual components 
of the photon in the transverse plane (the plane transverse to the collision axis) and their 
distribution in the photon can be obtained in the perturbative framework. The cross-section 
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a] ^(Q^, W"^) for the transversely {i = T) and longitudinally {i = L) polarised virtual photon 
is given by the following formula 

at%Q\W') = f^dz J d'r |vl>,(^, r)p O, (1) 

where r denotes the transverse separation between q and q in the colour dipole, z is the 
longitudinal momentum fraction of the quark in the photon and x is the Bjorken parameter, 
i.e. X = Q^/{2pq). The cross-section a{x,r^) is the (gg)dipoie^proton total cross-section, and 
|\E'j(z, r)p denotes the photon wave function squared and summed over the quark helicities, in 
the photon polarisation state indicated by i. The photon wave function is given by its quark 
flavour decomposition 

\nz,r)\' = Ei^f(^,^)r, (2) 



and 



I [z' + (1 - zf] e)Kl{efT) + m) Kl{efr), for i = T, 

= (3) 

' 4Q''z\l - zf K^itfr), for t = L, 



with 



ej = z{l-z)Q' + mj, (4) 

where Cf and rrif denote the charge and mass of the quark of flavour /. The functions Kq and 
Ki are the McDonald-Bessel functions. 

Equation (|l]) is in fact equivalent to the kt (or high energy) factorisation, which is the basic 
tool for calculating the observable quantities at low x B^, and the dipole-proton cross-section 



a{x,r ) is related to the unintegrated gluon distribution in the proton f(x, k ) ||33 



^(^, = ^ / - Ukrmx, k\ (5) 

In equation (^), k is the transverse momentum of the gluon and Jq{z) is the Bessel function. 
In the leading ln(l/x) approximation, the unintegrated gluon distribution f{x,k'^) is given by 
the solution of the BFKL equation which determines the 'hard' pomeron in the perturbative 
QCD ll^. The exchange of the perturbative QCD pomeron however violates unitarity at very 



large energies. The novel feature of the saturation model is the incorporation of the unitarity 
constraint on the level of the dipole-proton cross-section r^). This is achieved by imposing 
the saturation property, i.e. (t(x, r^) for r ^ Roix), where the saturation radius -Ro(^) is 

a decreasing function with decreasing x 

Rlix) ~ x\ (6) 

and the cross-section do is independent of x. In the limit r ^ 0, it follows from the perturba- 
tive QCD calculations, that the dipole cross-section exhibits the colour transparency property 
behaving as r^) ~ r'^/R^{x) (modulo a logarithmic correction which can modify the 
r dependence). Those two properties are economically summarised by the following simple 
parametrisation 

a{xy) = ao[l-exp(-rV(4i?^(x))]. (7) 
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^/\/\/\/\y<C^^^J^] 
1-Zi 




Figure 1: The diagram illustrating the 7*7* interaction in the dipole representation, see formula 



2.2 Generalisation of the GBW model for the two-photon cross- 
section 

The description of the 7*7* total cross-sections within the formahsm utihsing the interaction 
of two colour dipoles, which the virtual photons fluctuate into, has been discussed in detail 
The dipole-dipole cross-sections were assumed to be given by the Stochastic Vacuum 



m 



23-25 



Model p3|,^ or to follow from the BFKL formalism |2^. The novel feature of our approach 



is incorporation of the saturation property of the dipole-dipole cross-section. This makes it 
possible, in particular, to describe in a unified way the variation of the energy dependence of the 
cross-sections with the change of the virtualities of the photons. In terms of the virtual photon 
four-momenta qi and q2 we have Qi 2 = ~Qi2 = ill + ^2)^, see Fig. [|. The extension 

the saturation model to the case of Y iQi)!* (Q^) cross-sections for arbitrary virtualities QI2 
is given below. 

A formula for the two-photon cross-section part coming from the exchange of gluonic degrees 
of freedom reads 



1 1 

ajti^^ 7 Jo J ' 

(8) 

where the indices i,j label the polarisation states of the virtual photons, i.e. T or L. The 
wave functions ^"(2;fc,r) are given by equations (j^), with e'j defined by equation @ being 
replaced by (e^)^ = Zk{l — Zk)Ql + m^, k = 1,2 



and <Ja'i{xab,ri,r2) are the dipole-dipole 



total cross-sections corresponding to their different flavour content specified by the a and b. 

Inspired by the GBW simple choice for the dipole-proton cross-section, we use the following 
parametrisation of the dipole-dipole cross-section (Ja^b{xab,^i,f2) 



1 — exp 



cff 



^RoiXab) 



(9) 
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where for Xab we take the following expression symmetric in (1, 2) 



QI + QI + Ami + ^< 
W^ + Ql + Ql 



(10) 



which allows an extension of the model down to the limit Qi 2 = 0. Note, that Xab depends 
on the flavour of scattering quarks. We use the same parametrisation of the saturation radius 
Rq{x) as that in equation (7) in ||16|, i.e. 



A/2 



(11) 



and adopt the same set of parameters defining this quantity as those in ]T6|] . For the saturation 
value a^'^ of the dipole-dipole cross-section (cf. equation (|])) we set 



a,b 



(12) 



where cxo is the same as that in [0. For light flavours, equation (^) can be justified by the 
quark counting rule, as the ratio between the number of constituent quarks in a photon and 
the corresponding number of constituent quarks in the proton. Following |jl6|, we also use the 
same value of cxo'* for all flavours. 



Three scenarios for reff(ri,r2) are considered: 

+ ' 
min(r^, r|), 

3. = min(r^, r|)[l + ln(max(ri, r2)/ min(ri, r2))]. 

The first two cases are simple generalisations of the parametrisation adopted in the case of 
the dipole-hadron scattering, i.e. we get ~ r\ (r^g ~ r|) in the configurations r| ^ 
{rl ^ r|). Case (3) is motivated by the two-gluon exchange between the colour dipoles, giving 
the following cross-section 




[1 - Jo{kri)][l - Jo{kr2)] = min(r^, r2)[l + ln(max(ri, r2)/min(ri, r2))]. (13) 



In all three cases the dipole-dipole cross-section exhibits colour transparency, i.e. (jfl{x, ri, r2) 
for ri — > or r2 ^ 0. 



1. r, 



eff 



2. r'e. 



The formulae given above correspond to the 'pomeron' contribution to the 7*7* total cross- 
sections. This means that they represent the exchange of gluonic degrees of freedom giving rise 
to the component of cross-sections which, at high energies, does not decrease with increasing 
energy. Equation (H), defining the dipole-dipole cross-section, similarly to equation (^ for 
the dipole-proton cross-section, interpolates between the 'hard pomeron' at small transverse 
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separations ri 2 and the 'soft pomeron' at large transverse separations. In more detail, for 
small values of (i.e. for r^g ^ AE^{x)) one gets 

2 

<(x,,,ri,r2) ^ < \Zf- . ~ x-b. (14) 

which can be interpreted as the 'hard pomeron' contribution if the parameter A is identified 
with the 'hard pomeron' intercept. On the other hand for large dipole sizes we have 

f^aj(^afe,ri,r2) ~ (To'^ (15) 

i.e. the cross-section is only slowly varying with the energy in accordance with what is observed 
for the 'soft pomeron'. It should be emphasised that the structure of the saturation model is 
different from that corresponding to two separate (i.e. 'hard' and 'soft') pomeron contributions 



34| . The soft pomeron in the saturation model appears rather as an effect of unitarisation of 
the exchange amplitude of the 'hard' contribution, provided by multiple exchanges and self- 
interaction of the 'hard' pomerons. 



2.3 Scaling of the cross-sections 

The virtualities of the two photons can be arbitrary, thus the model can describe the following 
three cases of physical and phenomenological interest: 

1. The case Q\ = Q\ = ^ corresponding to the interaction of two real photons. 

2. The case Q\ ~ Q\ (with large Qi 2 ) corresponding to the interaction of two (highly) 
virtual photons. The relevant cross-section can be extracted from the measurement of 
double tagged events e^e~ —>■ e^e~ + hadrons. 

3. The case Q\ ^ Q\ corresponding to probing the structure of virtual {Ql > 0) or real 
{Q2 = 0) photon at small values of the Bjorken parameter x = Q\/ {2qiq2). For instance, 
the structure function F2 (x, Q'^) of the real photon {Q\ = 0,Ql = Q"^) is related in the 
following way to the 7*7 total cross-sections 

Fl{x, Q') = -^MW, Q\ = 0) + a^;r{W\ Q\ Q\ = 0)]. (16) 

47r ttem 

Let us now examine the high energy limit in all three cases and compare them with the 7*p 
case. We shall show that in the kinematical configurations in which the saturation effects are 
important, the two-photon cross-sections are more singular in the high energy limit than the 
7*p cross-section. To be precise, the high energy behaviour of the total 7*7* cross-sections will 
be enhanced by additional factors of the large logarithm \n[Q'^Rl{x)]. Also the behaviour of 
the 77 total cross-section will be shown to be enhanced by an additional power of ln(iy^/m^), 
i.e. the 77 total cross-section becomes steeper function of W than the 7p total cross-section, 
which seems to be confirmed experimentally. This enhancement will be a direct consequence 
of the singular behaviour of the photon wave function for small dipole sizes related to the 
point-like component of the photon. For simplicity, in the analysis given below, we focus on the 
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contributions to the two-photon cross-sections of the hght quarks m, d and s with the same mass 
rriu = md = = nig. The heavy flavour components exhibit the same general properties as the 
light flavour ones, but in the presently available kinematic range, the transition to saturation 
regime may not be observed. 

Let us first examine case (1) for the 77 total cross-section. The dominant contribution to 
the integrals in equation (|D comes from the region 4_Rq(x) <^ rf <^ r| -C l/'^^g, 4i?Q(x) <^ 
r2'^rf<^ where x = Sm^/W^ and Ro{x) is given by equation (pT]). In this region the 

short distance approximation of the (transverse) photon wave-function may be used 

\^t{z,v)\' ~ i, (17) 
and the corresponding contribution to the total 77 cross-section is 

a,,iW') ~ /'^'"' ^ ^ ~ \n'[ARl{x)ml] ~ In'iW'/W^). (18) 

This should be compared with the 7p total cross-section, where the saturation model extended 
down to the photoproduction limit gives a^p(W^) ~ \n{W'^ /Wq). 



Case (2) of the 7*7* cross-section in the configuration Ql ^ Ql ^ Q^, with Q"^ being large, is 
regarded as a very useful tool for probing the bare 'hard' pomeron exchange amplitude p5| , |36[ . 
The short-distance approximation (|l^ of the photon wave function is now valid in the region 

rlzk{l-Zk)Ql < 1, A; = 1,2. (19) 

The saturation model predicts different high energy behaviour of the 7*7* cross-section de- 
pending on whether > Q^{x) or < Q^x), with the saturation scale Q'^{x) ~ 1/Rl{x). 
Thus in the region > Ql{x) we get (modulo logarithmic corrections) 

a,,,, Q\-Q\Ql- Q') ~ yw^y (20) 

which exhibits the rise W'^^ and l/Q"^ dependence, characteristic for the hard pomeron ex- 
change. In the region < Q^{x) the 7*7* cross-section has the saturation property, i.e. 

a) the 1/Q^ behaviour is changed into a weakly varying function of Q^; 

b) the power-like W"^^ behaviour is replaced by a moderately increasing function of W"^. 

The leading behaviour of a^*^* for Q'^Rq{x) <^ 1 comes from the 'strongly ordered' configura- 
tions in integrals (^ defining the 7*7* total cross-sections 

> > > ARlix), (21) 



with k = 1 and 2. From equations (|8| JT7| , PT| ) one obtains 

a^*,*{W\Ql^Q\Ql^Q^) ~ \n\Q^Rl{x))[l + 0{l/ \n{Q^Rl{x))]. (22) 
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Note that in the saturation model at high energies, the 'pomeron' component of the 7*7* 
total cross-section is, to a good approximation, a function of only two variables ti = QfRKx) 
and T2 = Q2-^o(^)- analogous 'geometric scaling' was found in the DIS data A weak 



breaking of the scaling property in both cases occurs due to the presence of quark masses. For 
Qi^Ql^ Ti ~ r2 ~ r), the 7*7* cross-section exhibits the 1/r behaviour at large r 

(see eq. (|20|)) and reaches the saturation limit corresponding to a slowly varying function of r 
for small values of r. It should be observed that the leading behaviour at small r for the 7*7* to- 
tal cross-section (see equation (0)) is more singular than for the 7*p case, where a^.p ~ ln(l/r). 

Finally, in case (3) corresponding to probing the structure of the real (or quasi-real) photons 
at low X and large we find 

F^{x,Q') ~ x-\ (23) 

for Q'^ > Ql{x) and 

F^{x,Q') ~ Q'ln'iQ'Rlix)], (24) 

for g2 < g2 ^ Q2^^y 



2.4 Non-pomeron contributions 



In order to get a complete description of 7*7* interactions, which could be extended down to 
values of ly ~ 10 GeV, we should add to the 'pomeron' contribution defined by equation (H) 
the non-pomeron reggeon and QPM terms The additional contributions are characterised 
by a decreasing energy dependence, i.e. ~ for the reggeon and ~ (with InW 

corrections) for QPM. The QPM contribution, represented by the quark box diagrams, is 
well known and the cross-sections are given, for instance, in |3^. The reggeon contribution 
represents a non-perturbative phenomenon related to Regge trajectories of light mesons. It is 
known mainly from fits to total hadronic cross-sections and to the proton structure function F2. 
A state-of-art parametrisation of the reggeon exchange cross-section in two-photon interactions 
is given by the following expression 



a\W',QlQl) 



"em 



a2 



{a2 + Qi)ia2 + Qi] 



02 



(25) 



Originally, it was set [g^: A2 = 0.38, 02 = 0.3 GeV and rj = 0.45. However, those parameters 
were obtained with certain assumptions concerning the pomeron exchange which are different 



from ours. Therefore, it is legitimate to modify the parameters of |]2^ while retaining its 
functional form. We have in particular found that in order to get a good description of the 
data on the 77 total cross-section in the 'low' energy region, W < 10 GeV, one has to set 
T) = 0.3. This happens to be consistent with the recent observation that the intercept of the /2 
trajectory, which contributes to the two-photon cross- sect ion, can be expected to be equal to 
0.7 pOi. We fitted the other parameters, A2 and a2, to the data on two-photon collisions. 



Finally, note, that the decomposition of the reggeon term into different photon polarisation 
states has not been specified. In our analysis we assume, that the reggeon couples only to 
transverse photons. This arbitrary assumption, does not influence significantly our results for 
the studied observables. 
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2.5 Threshold corrections 



Strictly speaking, both the dipole model (accounting for the 'pomerons') and the Regge model 
of the total cross-section are formulated in the high energy limit x ~ Q'^/W^ <^ 1. When 
extending applicability of these models up to larger values of x, for instance x ~ 0.1, threshold 
correction factors should be taken into account. Namely, the cross-section should vanish when 
X — 1 as a power of 1 — x. In the case of 7*p scattering, the form of the cross-section at 
X ~ 1 is governed by the number of spectator quarks ngpect in the proton which do not interact 
directly with the photon. To be precise, it follows from the dimensional-counting rules that 
for a subprocess with a given number of spectators, at x ~ 1, the cross-section takes the form 
(T^*p(x, Q^) ~ (1 — xY^^p""^'^ (where the dependence is suppressed). A possible way to 
combine the small x dependence of the cross-section in the Regge model with the latter result 
is to include (1 — x)^"'=p'='''~^ as a multiplicative correction factor to the asymptotic cross-section 
from the pomeron or a subleading reggeon exchange For the pomeron exchange in 7*p 

scattering one has Aspect = 4 and for the other reggeons (/2 and 02), rigpect = 2. It is clear, that 
a similar procedure may be applied for the saturation model as well. 

In the case of two-photon collisions, one of the photons plays the role of the target, probed 
by the other photon. In the dipole representation, the number of valence quarks in the target 
photon equals two, to be compared with three valence quarks in the proton. Thus, for the non- 
pomeron reggeons one has ngpcct = 1 and for the dipole-dipole scattering component, which 
represents the 'pomeron' exchange, one obtains rigpect = 3. Recall, that in order to extend 
the saturation model to describe real photons, we use the variable Xab (see eq. (0)) instead 
of X. Also here, we represent the threshold correction factors using Xab- Thus we multiply the 
reggeon term (|25|) by (1 — Xgg) with Xgg obtained from eq. (Il0| ) with = = m^, that is 
defined by the light quark mass rUg. For the dipole-dipole scattering cross-section, Xab depends 
upon the flavour of quarks which span the dipoles. Hence, in our final formulae we multiply 
the dipole-dipole cross-section o"f'^(xafc, ri, r2), (see (^) by the factor (1 — Xa^)^. 



2.6 Final formulae 

For clarity, we collect the components of the two-photon cross-section presented above. The 
total 'y* {Qi)'y* {Q2) cross-section reads 

<*(W^^ Ql Ql) = ~4(W\ Ql Ql) + a^iW, Ql QI)6,tS,t + ^(1V^ Ql Ql), (26) 

where d'fj(W'^,Ql,Ql) is the gluonic component, corresponding to dipole-dipole scattering, as 
in eq. (|^), but with the dipole-dipole cross-section including the threshold correction factor 

afj(xa6,ri,r2) = (1 - x^;,) J(xa6, n, r2), (27) 

c.f. eq. (P), and Xab is given by eq. (jl^). The sub-leading reggeon contributes only to scattering 
of two transversely polarised photons and also contains a threshold correction 

a^{W^ Ql Ql) = (1 - x) a^iWl Ql Ql), (28) 

with 

_ _ QI + QI + 8m^ 

""'W^ + Ql + QV ^'^^ 
The third term a^J'^iW'^, Ql Ql) is the standard QPM contribution and is taken from |38 
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3 Comparison to experimental data 



3.1 Parameters of models 

In the comparison to the data we study three models, based on all cases for the effective 
radius, as described in Section 2.2. We will refer to the these models as Model 1, 2 and 
3, corresponding to the choice of the dipole-dipole cross-section. Let us recall that we take 
without any modification the parameters of the GBW model: ctq = 29.13 mb, Xq = 0.41 ■ 10~^ 
and A = 0.277. However, we fit the light quark mass to the two-photon data, since it is not very 
well constrained by the GBW fit, as we explicitly verified. On the other hand, the sensitivity 
of the choice of the mass appears to be large for the two-photon total cross-section. We find 
that the optimal values of the light quark (u, d and s) masses niq are 0.21, 0.23 and 0.30 GeV 
in Model 1, 2 and 3 correspondingly. Also, the masses of the charm and bottom quark are 
tuned within the range allowed by current measurements, to get the optimal global description 
in Model 1, r^g = rlrl/{rl + rl), which agrees best with data. For the charm quark we use 
rric = 1.3 GeV and for bottom rub = 4.5 GeV. Moreover, we re-fit t], A2 and 02 parameters in 
the reggeon term (^), which is legitimate because the 'pomeron' term which we use is different 
from this following from the model of two pomerons used in p4|. We find that the values 



r] = 0.3, A2 = 0.26 and 02 = 0.2 GeV^ give the best description of data, when combined 
with the saturation model. The values of masses listed above are consistently used also in the 
quark box contribution (QPM). The Models, which we shall mention from now on, contain the 
saturation models described in Section 2, combined with the reggeon and QPM contribution. 

It should be stressed, that most of the data relevant for our study were collected with the 



help of the comprehensive review [41 . 



3.2 The test case: the "fp total cross-section 

As stated above, we modify the quark mass of the GBW saturation model and the Donnachie- 
Dosch-Rueter parametrisation for the subleading reggeon. Certainly, one has to ensure this 
change does not spoil the quality of the GBW description. Besides that, it is necessary to 
check that the reggeon term with the modified exponent rj = 0.3 allows the extension of the 
GBW model for the jp total cross-section down to low values of ~ 3 GeV as well. Thus we 
calculated the dipole-proton scattering contribution using the original GBW approach, with 
the light quark mass, rriq, set to 0.21 GeV, as in Model 1, and added the reggeon term 

where A^p was fitted to data and the best value reads A^p = 0.135 mb. The result is given in 
Fig. 1^, where the cross-section from Model 1 is compared to the data, taken from Ref. and 



to the classical Donnachie-Landshoff fit [H2|. In the same figure we also show the decomposition 



of the total cross-section into the gluonic contribution, given by the saturation model and the 
reggeon component. Both contributions have been multiplied by a correction factor of the form 
(1 — a;)^"=p°'^*~\ as described in Sec. 2.5, with ngpect = 2 for the reggeon exchange, and ngpect = 4 
for the dipole-proton scattering. The fitted curve, with only one free parameter A^p follows the 
data accurately, suggesting that the model has certain universal properties. 
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1 10 100 

W [GeV] 

Figure 2: The total 7p cross-section - predictions from the GBW model with the light quark mass rriq 
set to 0.21 GeV and the charmed quark mass nic = 1-3 GeV, supplemented by the reggeon term 1^^), 
compared to data and to the Donnachie-Landshoff fit ^4^1 . Also shown are the gluonic and reggeon 
components of the full result in our model. The curves are cut atW = 3 GeV. 



3.3 Total 77 cross-section 

The available data for the 77 total cross-section range from the 77 energy W equal to about 



1 GeV up to about 160 GeV see Fig. ^ The experimental errors of the data are, 

unfortunately, rather large. One of the reasons is that those data were taken for virtual photons 
coming from electron beams and then the results were extrapolated to zero virtualities. Another 
problem appeared to be very important in LEP measurements where the incoming e"*" and 
e~, and a substantial fraction of the produced hadrons go into the beam pipe and cannot be 
detected. Extraction of the actual 77 collision energy is therefore needed from the visible energy, 
which is a model dependent procedure and introduces large systematic errors. In particular, 
it is well known, that the data for the 77 total cross-section from LEP depend on the Monte 
Carlo method applied for the unfolding. 

The data interpreted with Pythia tend to be larger and exhibit a steeper rise with W 
than those unfolded with P HO JET (see Fig. ^). In the saturation model it is difficult to 
obtain a cross-section consistent with the Pythi A- unfolded data. Besides that, Phojet is a 
Monte Carlo program dedicated to describe two-photon interactions and the description of the 
crucial hadron emissions close to two-photon collision axis, is elaborated in more detail. Thus 
we choose to follow the Phojet unfolded data in our analysis. In a more conservative approach 
one should include the difference between the cross-sections unfolded with different programs 
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Figure 3: The total 77 cross-section: (a) data confronted with predictions from all three Models and 
(h) contributions to the Model 1 result: gluonic, reggeon and QPM. 



into the systematic error. This would only make the data less constraining and would not spoil 
the quality of the fit. 

In Fig. ^ we show the total 77 cross-section from the Models, obtained using eq. ( p6|) with 
i = j = T. The agreement with data is very good down to ~ 3 GeV. It is interesting 
to observe, that the Models strongly favour the Phojet unfolded data, and that the energy 
dependence of the total 77 cross-section (PHOJET-unfolded) at high W is very well reproduced 
by all three Models, see Fig. 03. Recall, that the steeper VT-dependence found in the two- 
photon cross-section, as compared to the hadronic and the photoproduction cross-sections, is 
naturally explained by the presence of additional factors of In^ W and In W respectively, as 
discussed in Section 2.3. Predictions for W in the range to be probed in future linear colliders 
(i.e. W < 1 TeV) are stable against variations of the details of the saturation models, provided 
that the models are adjusted to fit the available data. 
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Figure 4: Tota/ 7*7* cross-section /or (a) = 3.5 GeV^ (b) = 14 GeV^ and (c) = 17.9 QeV^ 
- comparison between LEP data and the Models plotted as a function ofY = \n.{W'^ / Q"^) . Also shown 
is the result of Ref. Jsdjl based on the BFKL formalism with subleading corrections, supplemented by 
the QPM term, the soft pomeron and the subleading reggeon contributions. 



3.4 Total 7*7* cross-section 



The data |^8|,^ for the total 7*7* cross-section are extracted from so-called double-tagged 
events, that is from e^e~ events in which both the scattered electrons are measured and hadrons 
are produced. In such events measurement of the kinematical variables of the leptons determines 
both the virtualities and Q2 of the colliding photons and the collision energy W. The tagging 
angles in LEP experiments restrict the virtualities to be similar, i.e Qf ^ Q2 = Q^- The data 
are available from LEP for average values = 3.5 GeV^, 14 GeV^ and = 17.9 GeV^ in a 
wide range of W. 

In Fig. ^,b,c those data are compared with the curves from the Models. As an estimate 
of the total 7*7* cross-section we use a simple sum of the cross-sections cr-°* (eq. (p6D) over 
transverse and longitudinal polarisations i and j of both photons. In addition we plot also 
the prediction obtained in Ref. |^ by solving the BFKL equation with non-leading effects, 



and added phenomenological soft pomeron and reggeon contributions and the QPM term. The 
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Figure 5: The total cross-section - the decomposition of Model 1 into QPM, gluonic and reggeon 
components for (a) low virtuality = 3.5 GeV^ and (h) large virtuality = 17.9 GeV^. 



latter prediction was found to describe the measured e~^e~ differential cross-section for hadron 
production in double tagged events . As can be seen, Model 1 fits data as well as the result 
based on the BFKL solution. Model 2 is slightly worse than Model 1 and Model 3 may be 
rejected. 

Since the virtuality is high, the unitarity corrections are not important here. For the 
same reason, the results are not sensitive to the choice of the quark masses and the parameters 
of the reggeon term. As seen in Fig. ^, where the components of the 7*7* total cross-section 
from Model 1 are plotted, the cross-section is dominated by the QPM and the 'pomeron' 
contributions. Moreover, the perturbative approximation for the photon wave function is fully 
justified in this case. Thus, in this measurement the form of the dipole-dipole cross-section is 
directly probed. 



3.5 Photon structure 

The quasi-real photon structure may be probed in single tagged e+e" events. In this case 
one of the electrons scatters with a larger momentum transfer Qf = which corresponds 
to the emitted photon virtuality and the other electron scatters at a low angle, producing 
predominantly a virtual photon with very low virtuality Q| ~ 0. Thus, the measurement of the 
cross-section for the 7*(<5i)7*(<52) collision at the energy W can be used to extract the almost 
real photon structure function (x, Q^), with x = Q"^/ (W^ + Q^), see eq. (|16D. 

Essentially, the parameters of the models are constrained by the data for the total 77 and 
7*7* cross-sections so here we are presenting a parameter free result. In Fig. ^ we show the 
comparison of our predictions with the experimental data |^0|,^ for the virtuality in the 



range from (a) 1.9 to 2.8 GeV^, (b) 3.7 to 5.1 GeV^, (c) 8.9 to 12.0 GeV^ and finally (d) from 
16.0 to 23.1 GeV^. Note, that in each plot the data for various virtualities are combined, 
which may give rise to systematic effects, see for instance Fig. ^3. In each plot the value of 
virtuality adopted to obtain the theoretical curve is indicated and was selected to match 
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Figure 6: The photon structure function F2{x,Q'^): the experimental data compared to predictions 
following from the Models for various : (a) from 1.9 to 2.8 GeV^ , (h) from 3.7 to 5.1 GeV^ , (c) from 
8.9 to 12.0 GeV^ and (d) from 16.0 to 23.1 GeV^ . The curves are cut at values of x corresponding to 
W = 2, GeV. 



the average value of the data-set containing the best data at low x. We also show in Fig. ^ 
the contributions to F2 in Model 1. As already stated, the importance of the reggeon term is 
not very large and decreases with increasing Q^. 

Model 1, favoured by the 7*7* data provides the best description of F2 as well. In the 
region of a; > 0.1 the agreement of the Model 1 with data is surprisingly good which was not 
a priori expected from the model based on the large energy approximation. Note however, that 
a systematic tendency occurs for all the Models to overestimate the data for larger Q^. 

It is straightforward to obtain in this framework predictions for the virtual photon (with 
Q\ = P^) structure function F2 (x, Q^; P^) in the low x domain. However, there exist only very 
few data on this observable so we do not present our predictions for this quantity. Nevertheless, 
possible experimental study of F2 (x, Q^;P^) would certainly provide another interesting test 
of the saturation model. 
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a) b) 

Figure 7: The photon structure function F2{x,Q'^) from Model 1 compared to the experimental data 
for (a) low 1.9 to 2.8 GeV^, and (h) large from 16.0 to 23.1 GeV^. 



3.6 Heavy flavour production 

Another interesting process which we have studied in the dipole model is the production of 
heavy flavours (charm and bottom) in 77 collisions. Heavy quarks can be produced by three 
mechanisms: 



1. The direct production in which both photons couple to the same heavy quark line, which 
corresponds to a component of the quark box diagram (see Fig. ^). 

2. The direct photoproduction off the resolved photon, which would involve a fluctuation 
of one of the photons into a heavy quark-antiquark pair and scattering of the pair off 
the other, resolved photon by exchange of gluons. This phenomenon is accounted for 
in the dipole model and the cross-section may be obtained by restricting the sum over 
the flavours in eq. (0) to the case, in which at least one dipole is composed of the heavy 
quarks only (see Fig. ^). 

3. The hard fragmentation and rescattering contributions. The first one corresponds to 
production of a heavy quark pair in the fragmentation process of an initial light quark 
pair. The initial pairs can be produced either through the box diagram (Fig. ^) or as the 
colour dipoles (Fig. ^) which are present in the model. Note also that in the saturation 
model, in the case of real photons and the original dipoles composed of light quarks, 
abundant rescattering of cascading gluons occur in which heavy quarks may be produced. 
This would be the other, rescattering mechanism. The estimate of such effects is, so far, 
beyond the reach of our model and we do not take into account these contributions. 



The reggeon exchange is a non-perturbative phenomenon and should not contribute to heavy 
fiavour production, so it is assumed to vanish here. In Fig. ^ we plot the predictions from all 
three Models compared with L3 data on charm production ||5^. The best model. Model 1, is 
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Figure 8: Diagrams illustrating contributions to the heavy quark production, at the amplitude level: 
(a) heavy quark box diagram, (b) direct production of heavy flavours in one of the dipoles, (c) example 
of production of a heavy quark pair through a hard fragmentation process in a box diagram, and (d) 
representation of production by hard fragmentation from a light quark dipole, or by gluon rescattering. 



below the data leaving some room for a possible contribution from the fragmentation. The 
shape of the cross-section is well reproduced. 

Production of bottom quarks in two almost real photon collisions was investigated ex- 
perimentally by the L3 |^| the OPAL |54] collaborations. There, the measured process 



was e'^e e^e bbX, with anti-tagged electrons at e~^e invariant collision energies -y/Sgg 
between 189 GeV and 202 GeV. The total cross-section for this reaction was found to be 
13.1 ± 2.0 (stat) ± 2.4 (syst) pb (L3) and 14.2 ± 2.5 (stat) ± 5 (syst) pb (OPAL) whereas the 
theoretical estimate from Model 1 for \/s^^ = 200 GeV gives about 5.5 pb with less than 10% 
uncertainty related to the choice of 6-quark mass. This is significantly below the experimental 
data but above the expectations of 3 ± 1 pb (see e.g. [@]), based on standard QCD calculations 
with the use of the resolved photon approximation. 

In conclusion, the saturation model underestimates the cross-section for production of heavy 
quarks and the discrepancy increases with increasing quark mass, or perhaps, decreasing electric 
charge. This may be a hint that the fragmentation and rescattering mechanisms of heavy quark 
production are, indeed, important. 
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Figure 9: The cross section for the inclusive charm production in 77 collisions: (a) results for all 
three Models and (b ) the decomposition of the result from Model 1 on the QPM and gluonic component. 



4 Predictions for future colliders 

Two-photon processes will be important at possible future e~^e~ colliders, like TESLA, where the 
available photon-photon collision energy might reach 500 GeV or even 1 TeV |^ . Thus, we give 
predictions from Model 1 for the energy dependence of 7*(<5^)7*(<5^) , the gluonic component 
a^.^,{W^,Q'^,Q'^) (defined by eqs. (D and (0)) of the total cross-section a'^'l^^QlQlW) 
for Ql = Ql = Q^. In Fig. p!o|a,b we show the results in terms of a re-scaled quan- 
tity Q'^a^,^,{W'^,Q'^,Q'^) , with = max(Q^4m2) for various between and 10 GeV^ 
(Fig. p!0|a) and between 10 and 200 GeV^ (Fig. |lO|b). Of course, the cross-section for gluon 
exchange cr^,^*(W^^, Q^, Q"^) has to be combined with the standard QPM and reggeon terms in 
order to get a complete description of the total cross- sect ion, as described in Sec. 2.6. 

For comparison the same quantity, but obtained from the solution of the BFKL equation 
with subleading corrections in the perturbative domain is presented in Fig. ^c. Note, 
that the latter result contains not only the perturbative BFKL part but also the soft pomeron 
contribution, obtained using the Regge factorisation. Recall that the relative importance of the 
soft pomeron term quickly decreases with the increasing photon virtuality Q^. 

It is interesting to observe, that in the saturation model, at higher W where the threshold 
corrections are negligible, the quantity Q'^a^,^*(W'^,Q'^,Q'^) increases with Q"^ up to about 
10 GeV^, see Fig. |rU|a and for > 10 GeV^ it decreases with Q^. The reason for the increase for 
smaller is the rising contribution from the charmed quark dipoles. The relative suppression 
of the charmed quark contribution at = 0, in comparison to the light quarks due to the 
higher charm mass, is becoming less important towards higher > Arn^, when the typical 
scales in those two cases become similar. 

On the other hand, the bottom quark has a relatively small charge of 1/3 e and such thresh- 
old effects are much less pronounced. Thus, for Q"^ > 4m^ the cross-section should enter the 
geometric scaling regime. The unitarity corrections may be neglected at higher and it fol- 
lows that Q^CTy,^. (W^^, Q^, Q^) ~ {W^/Q"^)^, modulo threshold corrections. This is the reason 
why one sees a monotonical decrease of Q'^d-^,^,(W'^, Q^, Q"^) with Q"^ in Fig. p!0|b. 
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Figure 10: Predictions for the re-scaled gluonic component of the two virtual photon cross-section 
Q'^a^t^t{W'^,Q'^,Q'^) obtained using the parameters of Model 1 for (a) between and 10 GeV^ and 
(b) between 10 and 200 GeV^ , and (c) the same quantity from the BFKL solution combined with the 
soft pomeron contribution. 



For Q in the perturbative domain, the difference between the results from saturation model 
and the BFKL predictions is not large but grows with Q"^. The tendency of Q'^a^*^, (W^, Q^, Q^) 
to decrease with increasing in the BFKL approach may be traced back to important thresh- 



old effects and the running of the QCD coupling, which were taken into account in |3B . 

Let us also recall that the prediction for the 77 total cross-section for W in the TeV range is 
not sensitive to the choice of the form of the dipole-dipole cross-section, see Fig. ^. However, it 
does rely on the accuracy of the data unfolding with Phojet. Thus, the systematic uncertainty 
of the data unfolding at LEP propagates into the model predictions. 

In Fig. |rT| the real photon structure F2 (x, Q^) is plotted for various between 1 GeV^ 
and 200 GeV^ for x down to 10~^. The kinematical range was chosen to be relevant for the 
future linear collider measurements. For completeness, in Fig. we also give the dependence 
of the cross-section for heavy quarks production in 77 collisions in a wide range of two-photon 
collision energy W. We indicate the effect of the quark mass variation both for charm and for 
bottom. One should, however, keep in mind, that the saturation model in the present form 
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Figure 11: The real photon structure function F2{x,Q'^) for various from Model 1. 



gives slightly too low cross-sections for c quarks production and significantly too low (by factor 
of about 2-2.5) for b quarks. 

5 Conclusions 

In this paper we have extended the saturation model proposed by Golec-Biernat and Wiisthoff 
for 7*p cross-sections, to describe two-photon processes at high energies. This extension re- 
quired an explicit model for the scattering of two colour dipoles. We considered three models 
of this cross-section, all of them exhibiting the essential feature of colour transparency for small 
dipoles, and the saturation property for large ones. We kept the GBW form of the unitarising 
function and the original parameters, except for changing the values of quark masses, which 
was necessary to describe the data on the total two real photon cross-section. We have ex- 
plicitly checked that these modifications had not spoiled the fit to photoproduction data. In 
order to obtain a more complete description applicable at lower energies the saturation model 
has been combined with other, well known contributions related to the quark box diagram and 
non-pomeron reggeon exchange. Those mechanisms become dominant when the collision en- 
ergy is comparable to photon virtualities and quark masses. Standard multiplicative threshold 
correction factors, relevant for lower energies, have been included into the saturation model 
and the sub-leading reggeon contribution. 

Wc have analysed general features of the saturation model for the 77 scattering. We ob- 
served, that such models gave the energy (W) dependence for two-photon total cross-section in 
the saturation regime steeper than the photon-hadron and hadron-hadron total cross-sections, 
by additional factors of In W and In^ W respectively. In the non-saturated regime, a typical 
power-like growth W'^^, A ~ 0.3, with energy was obtained, characteristic for the hard pomeron 
exchange. 
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Figure 12: Cross-sections for heavy quark production in two real photon collisions from Model 1. 



The results from the studied models were compared with the data for different two-photon 
processes at high rapidity values: the total 77 cross-section, the total 7*7* cross-section for 
similar virtualities of the photons, the real photon structure function F2 and heavy flavour 
production. Free parameters of the models were fitted to the data. It was found, that the data 
favour one of the models for dipole-dipole cross-section, namely Model 1 presented in Section 
2.2 and Section 3.1. With this model a reasonable global description of the available two- photon 
data was obtained, except for the 6-quark production. Predictions for energies accessible at 
future linear colliders were formulated. It is however encouraging, that Model 2 (see Sec. 2.2 
and Sec. 3.1), being a significantly different generalisation of the original GBW model, gives 
the results close to Model 1, with a relative difference of less than 15%. This means that the 
sensitivity of the predictions from the saturation model to the details is not very significant. 

In summary, the saturation model was found to provide a simple and efficient framework 
to calculate observables in two-photon processes. This success supports strongly the idea of 
rapidity dependent saturation scale and improves the understanding of two-photon physics. 
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